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ABSTRACT 


We  have  studied  the  group  and  phase  speed  structure  of  long-range  acoustic  arrivals  using  linear  acoustics.  Here 
we  present  a  preliminary  study  of  some  of  these  same  paths  but  now  we  consider  a  higher  amplitude  source 
necessitating  a  nonlinear  acoustic  treatment.  This  research  involves  the  eventual  coupling  of  the  Nonlinear 
Progressive  Wave  Equation  (NPE)  with  an  Adiabatic  Normal  Mode  (ANM)  Model.  We  use  the  ANM  model 
(which  has  limited  range-dependent  capability)  because  our  previous  treatment  for  the  linear  problem  was  very 
revealing  with  respect  to  the  physics  of  the  arrival  structure. 
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OBJECTIVE 


The  objective  of  this  work  is  to  derive  new  technologies,  based  on  the  physics  of  underwater  sound  propagation 
(both  linear  and  nonlinear),  for  localizing  underwater  sound  sources  using  data  from  single  hydrophone  stations 
that  are  part  of  the  International  Monitoring  System. 
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Introduction 


According  to  linear  acoustic  theory,  the  arrival  structure  of  a  pulse  propagating  a  long  distance  in  the  ocean  is  a 
result  of  the  relationship  between  group  speed  and  phase  speed  of  its  various  components,  the  latter  often 
characterized  as  the  frequency  dependent  normal  modes,  or  more  simply,  acoustic  multipaths  (Kuperman  et  al, 
2001).  In  particular,  it  is  the  appropriately  integrated  group  speed  over  the  individual  acoustic/geographic  paths 
that  are  important,  and  the  arrival  structure  therefore  contains  source  property  information  including  its  location. 
This  group-vs-phase  speed  property  is,  in  linear  acoustics,  independent  of  source  amplitude.  Even  though  the 
nonlinear  or  amplitude-dependent  regime  becomes  more  complicated,  the  original  nonlinear  nature  of  the 
source  may  remain  as  a  scar  in  the  received  long  distance  linear  arrival.  The  goal  of  this  research  is  to  address 
the  possibility  of  this  potential  diagnostic.  This  paper  outlines  some  of  the  basic  linear-vs-nonlinear  issues  and 
presents  some  initial  calculations  illustrating  these  points. 

Group  Speed  vs  Phase  Speed  in  Linear  Acoustics 

Though  global  acoustic  propagation  can  occur  only  in  deep  water,  we  must  include  coastal  or  shallow  water, 
bottom  reflecting  environments  because  of  the  potential  location  of  a  high  amplitude  source.  The  classic 
group/phase  speed  representation  in  shallow  water  is  the  so-called  Pekeris  curves  (Jensen  et  al,  1994)  shown  in 
Figure  1  for  modal  propagation: 

We  are  basically  interested  in  the  regime  to  the  right  of  the  minimum  (Airy  phase)  of  the  group  speed  curves. 
Here  we  see  that  when  phase  speed  increases,  group  speed  decreases  corresponding  to  a  positive  waveguide 
invariant  fi  (Chuprov,  1982;  Brekhovskiskh  and  Lysanov,  1991;  D’Spain  and  Kuperman,  1999).  Physically,  this 
can  be  easily  explained  in  terms  of  paths  with  higher  angles  with  respect  to  the  horizontal,  and  therefore,  large 
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(horizontal)  phase  speeds,  propagate  energy  at  lower  speeds. 

On  the  other  hand,  the  set  of  deep  water  paths  shown  in  Figure  2  have  more  diversity  in  their  group-vs-phase 
speed  dependence.  Bottom  bounce  paths  have  the  same  relation  as  above,  but  the  deep  refracting  paths,  such  as 
the  RR  and  RSR  paths  typically  have  the  opposite  group-vs-phase  speed  dependence.  That  is,  group  speed 
increases  with  increasing  phase  speed  (|3  <  0)  because  the  deep  paths  refract  in  water  column  regions  of  higher 
speed  as  opposed  to  the  bottom  reflection  case  in  which  a  reflection  occurs  before  a  path  can  reach  a  sufficiently 
higher  speed  region  in  which  the  net  effect  over  the  path  is  to  increase  the  group  speed.  In  realistic  ocean 
environments,  the  sound  speed  varies  with  range  such  as  the  climatology  data  associated  with  recent  ATOC 
experiments  (Worcester  et  al,  1991;  ATOC  Consortium,  1998;  Worcester  et  al,  1999;  Colosi,  et  al,  1999; 

Heaney  and  Kuperman,  1998)  shown  in  Figure  3.  We  must  then  calculate  the  effective  group  speed  over  the 
whole  path.  An  adiabatic  normal  mode  computation  (Jensen  et  al,  Heaney  et  al,  1991)  is  often  adequate,  where, 
in  addition,  we  approximate  the  phase  speed  to  be  a  continuous  variable  because  of  the  typically  large  number 
of  modes,  whereas  in  reality  there  is  a  discrete  distribution  of  modes. 


The  local  phase  speed  of  the  nth  mode  at  the  source  location  (which  corresponds  to  a  launch  angle)  is  related  to 
the  modal  wavenumber  k„, 

Cp„(®)  =  f,  (1) 

where  to  is  the  angular  frequency  of  the  propagating  acoustic  mode.  In  a  range-dependent  environment,  we 
must  derive  an  effective  average  modal  group  speed,  u„,  from  the  range-averaged  slowness.  That  is,  we  note 
that  the  pulse  observables  at  range  r  are  arrival  times  that  are  given  by  r/u„.  Hence,  for  a  range-dependent 
environment  we  are  concerned  with  l/u„=  s  g„  where  the  overbar  denotes  range  averaging  and  Sg„  is  referred  to 
as  the  (local)  modal  group  slowness.  Since  the  adiabatic  mode  approximation  has  a  range-averaged  modal 
wavenumber  associated  with  each  term,  it  follows  that  (Jensen  et  al,  1994) 


Sg„ir,o}) 


(2) 


and  the  range-dependent  effective  group  speed  is, 

1 

<u„ir,co)>=Sg„{r,aj)  (3) 


Hence,  the  range-dependent  effective  group  speed  is  obtained  from  the  harmonic  average: 


1  1  V'  Ar' 

- 1 - —  (4) 

{u„(r,oj))  r^M„(r',c«) 

where  the  sum  is  taken  over  the  approximate  n  range-independent  subintervals,  Ar  « ,  of  the  propagation  path. 
Finally,  we  note  that  we  can  interchange  the  partial  with  the  integral  sign  in  Eq.  2  and  can  use  the  modal 
formula  for  the  (local  or  range-independent)  modal  group  speed  (Jensen  et  al,  1994), 


dk„  1  m  p^n{a>',z) 

dm  u„(a))  k„(m)'’o  p(^z)c^(z) 


(5) 


where  />„(co;  z)  are  the  normal  mode  eigenfunctions  of  the  pressure  field  in  the  ocean  waveguide  and  p(z)  and 
c(z)  are  the  density  and  sound  speed  profiles  as  a  function  of  depth. 

Applying  Eq.  (4)  to  the  data  in  Figure  3,  we  obtain  the  effective  group  speed  (at  a  final  specified  range)  vs 
phase  speed  (g-v)  curves  parameterized  by  frequency  shown  in  Figure  4.  These  curves  are  actually  a  summary 
of  the  expected  pulse  arrival  structure.  Thus,  for  example,  we  should  see  three  dominant  features  in  the  arrival 
structure.  At  the  lowest  phase  speed,  we  note  that  all  the  modal  group  speeds  are  independent  of  frequency  and 
these  components  have  the  same  and  lowest  group  speed.  Hence,  there  should  be  a  high-intensity  last  arrival 
since  all  components  arrive  at  the  same  time.  Next  there  should  another  high-intensity  arrival  (RR/RSR)  region 
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corresponding  to  a  phase  and  group  speed  of  (1516  km/s,  1483.6  km/s).  At  this  coordinate  there  is  a  stationary 
region  for  a  group  of  modes  (Kuperman  et  al,  2001).  Finally,  there  is  the  region  on  the  upper  right  that  cuts  off 
(to  dispersive,  lossy  bottom  interacting  modes.  This  high  group  speed  region,  which  describes  the  fastest  paths 
(earliest  arrivals)  should  be  low  in  amplitude.  The  above  is  descriptive  of  the  experimental  results  for  this 
ATOC  path  shown  in  Figure  5. 

Finally,  with  respect  to  linear  acoustics,  we  mention  that,  as  described  above,  since  shallow  and  deep  water 
propagation  have  different  g-v  curves,  a  transition  from  shallow  to  deep  water  will  have  a  profound  effect  on 
the  arrival  structure.  This  can  be  seen  from  Eq.  (4),  which  states  that  the  effective  group  speed  is  an  average 
over  group  slowness.  That  is,  group  speeds  add  like  parallel  resistors  and  so  a  small  continental  shelf  region 
will  have  a  large  effect  on  the  arrival  structure.  Hence,  we  have  a  linear  diagnostic  for  sources  originating  in 
shallow  vs  deep  water  regions. 

Nonlinear  Propagation  Effects 

Large  amplitude  sources  produce  nonlinear  effects  in  which  super-position  principle  breaks  down  and  new 
frequencies  are  generated.  These  effects,  even  after  decay  to  the  linear  acoustic  regime,  remain  as  possible 
diagnostics  of  high-energy  sources.  In  this  section,  we  summarize  some  of  the  nonlinear  waveguide  effects  that 
we  expect  will  be  useful.  We  will  base  the  discussion  on  the  Nonlinear  Progressive  Wave  Equation  (NPE) 
(McDonald  and  Kuperman,  1985;  McDonald  and  Kuperman,  1987). 

The  NPE  is  based  on  a  first-order  nonlinear  wave  equation  bearing  the  same  relation  to  hydrodynamic 
compressional  waves  that  the  Korteweg  -  de  Vries  (KdV)  equation  bears  to  weakly  nonlinear  incompressible 
waves  in  water.  The  NPE  is  cast  in  a  pulse-tracking  frame  of  reference  translating  at  a  speed  approximating  the 
speed  of  sound  (see  Figure  6),  and  includes  lowest  order  compressional  nonlinearity  along  with  diffraction  and 
refraction.  For  appropriate  environments,  it  possesses  similarity  solutions  (McDonald  and  Ambrosiano,  1988) 
analogous  to  the  solitary  wave  solutions  of  the  KdV  equation. 

Normally  we  use  pressure  in  linear  acoustics,  but  it  has  been  convenient  to  derive  the  NPE  in  terms  of  density; 
let  7?/(r,  t)  be  a  linear  baseline  acoustic  solution,  then  the  NPE  for  quadratic  nonlinearity  is  given  by 

(^0  ~  ^1  ’  (^) 

2r 

where,  oc  is  the  coeffcient  of  nonlinearity.  With  the  modal  sum  expressed  in  the  time  domain 

R\  {r,z,  t)  = '^  A „  (z)  cos(k„r  -mt  +  (p„),  (7) 

n 

we  obtain,  after  substantial  mathematical  manipulation,  a  nonlinear  modal  propagation  equation  for  the  depth 
integrated  field, 

/  dzR^=a^  Q(d„)sm2(k„r'+  (p„)  (8) 

n 

where  the  actual  form  of  Q  is  given  elsewhere  (McDonald,  2002).  The  important  point  here  is  the  generation  of 
harmonics.  This  result  is  a  particular  simple  form  obtained  by  restricting  the  background  environment  to  an 
analytical,  convergence  zone  generating  sound  speed  profile.  Here  we  begin  to  see  solutions  that  exhibit  the 
modal  behavior  that  must  accompany  nonlinear  acoustic  propagation;  that  is,  the  emergence  of  harmonics.  This 
idealized  result,  shown  in  Figure  7,  indicates  that  the  vertically  integrated  second  harmonic  maximizes  near  odd 
number  convergence  zones  and  represents  an  example  of  the  expected  departure  from  linear  results  of  the 
previous  section. 

The  nonlinear  tendency  toward  greater  bottom  interaction  has  been  documented  elsewhere  (Ambrosiano  et  al, 
1990)  using  the  NPE.  This  would  alter  the  modal  excitation  spectrum  toward  bottom  interacting  acoustic 
modes.  That  is,  wavefront  steepening  shifts  the  source  spectrum  toward  higher  frequency,  while  Mach 
reflection  increases  bottom  interaction.  A  schematic  modal  frequency  -  arrival  time  diagram,  shown  in  Figure 
8,  illustrates  that  bottom  interaction  gives  rise  to  clustered  low-frequency  arrivals  before  the  final  deep  water 
axial  cutoff  (the  lowest  group  speed).  Figure.  9  shows  a  plot  for  another  ATOC  example  to  be  considered  in  the 
context  of  Figure  8.  Thus,  one  evidence  of  a  nonlinear  event  in  the  water  column  could  be  the  appearance  of 
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another  focal  spot  at  low  frequency  in  Figure  9  corresponding  to  the  arrival  of  the  bottom  interacting  modes  vs 
the  already  existing  surface  interacting  focal  region  at  2010  seconds.  This  suggests  that  we  may  be  able  to  fold 
nonlinear  effects  into  a  g-v,  albeit  amplitude  dependent,  type  summary  of  the  propagation  conditions. 


CONCLUSIONS  AND  RECOMMENDATIONS 


We  have  shown  that  the  group  speed  dependence  on  phase  speed  (g-v  curves)  is  diagnostic  of  the  arrival 
structure  of  a  long  -distance  propagating  pulse  governed  by  linear  acoustics.  Further,  we  have  shown  that  the 
effective  group  speed  results  from  averaging  slowness.  The  different  g-v  structure  in  shallow  and  deep  water 
will  show  up  in  the  received  arrival  structure  of  a  pulse  traveling  between  shallow  and  deep  water.  Finally,  we 
have  also  demonstrated  that  nonlinear  modal  propagation  redistributes  the  modal  acoustic  energy  and  hence  the 
integrated  group  speeds.  Our  next  goal  is  to  cast  these  nonlinear  results  in  the  analogous  g-v  structure  we  have 
set  up  for  the  linear  case. 
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Figure  2:  Schematic  of  deep  water  paths. 
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Figure  3:  Sound  speed  profiles  for  a  3500  km  ATOC  propagation  path  as  obtained  from  a  climatology  data 
base. 
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Figure  4:  Group  speed  vs  phase  speed  (g-v)  curves  for  the  climatology  derived  sound  speeds  of  Figure  3. 


Figure  5:  ATOC  data:  The  pulse  arrival  structure  corresponding  to  Figs.  4  and  5. 
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Figure  6:  The  geometry  for  the  Nonlinear  Progressive  Wave  Equation  (NPE) 


Verttcally  Integrated  2nd  Harmonic  Q(k) 


Figure  7:  Second  harmonic  of  a  high  amplitude  source  for  the  particular  case  of  a  parabolic  slowness  profile 
has  local  maxima  near  odd  number  convergence  zones 


41 


23^  Hiif.cuS'ch  ftfi  fi'-n'r  Jftfj'tiinilit;  Mr'/iiit-fVtnff  t>J  iVdffA'rti^  ExpU>-i,t(tTi\  -  Ontiifbiit  2-S,  2^tf 


fh.rili.ini  |iik™r<ri>¥  [K^pirhiirii 


Figure  8:  Dispersion  curves  for  a  set  of  bottom  interacting  modes 
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Figure  9:  Modal  dispersion  and  arrival  structure  for  an  ATOC  propagation  example. 
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